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Abstract 

We establish exact results for the one-dimensional spin-orbital model for special values 
of the four-spin interaction V and dimerization parameter 5. The first exact result is at 
5 = 1/2 and V = —2. Because we find a very small but finite gap in this dimerized 
chain, this can serve as a very strong test case for numerical and approximate analytical 
techniques. The second result is for the homogeneous chain with V = —4 and gives 
evidence that the system has a spontaneously dimerized ground state. Numerical results 
indicate that the interplay between dimerization and interaction could result in gapless 
phases in the regime < y < — 2. 
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The renewed interest in spin-orbital models ^ Q stems from the belief that they are relevant 
to explain unusual magnetic properties of a wide range of materials. Representative examples 
are quasi-one-dimensional TDAE-Ceo material the Na2Ti2Sb02 and NaV205 compounds 

§], and two-dimensional LINIO2 samples [Q. In this letter we study the dimerized spin- 
orbital chain with Hamiltonian 

i=l = 1,2 2 = 1 

where a = 1,2 are two commuting sets of s = 1/2 operators acting on site i of a lattice 
of size L. Notice that the dimerization S is completely staggered, which seems to be the most 
interesting case at least in the context of spin ladders |[. The the four-spin coupling V in 
spin-orbital models originates in the standard superexchange interaction |]T|. In spin ladders, 
it can be interpreted as a result of the Coulomb interaction between holes in a doped phase 
0. Here, unless stated otherwise, we assume an antiferromagnetic J > Heisenberg energy 
scale. 

Much of the theoretical efforts have been focused on the study of the homogeneous chain 
both numerically [Tl|, 12| and analytically [13, 14]. In general, the interesting features occur 



in the strong coupling regime where perturbation theory does not work. In these cases it is 
important to have access to exact results since they play the role of relevant testing devices 
for approximate non-perturbative methods. One such example is the point V/J = 4m the 
homogeneous chain, which is equivalent to the integrable SU(4) exchange spin chain |T^. This 
fact has been not only useful to check numerical analysis |^, [l^ but also relevant to trigger 



further non-perturbative studies [|T5|. One of the purposes of this paper is to point out that 
there exists two other points in which the model (1) is still exactly solvable. Equally important, 
both results are concerned with the strong coupling regime \V\ » J, namely 

(I) 5 = 1/2, V/J = -2 and (II) S = 0, V/J = -4 (2) 



The exact solution enables us to derive analytical expressions for the ground state energy 
and the low-lying gap excitation. To our knowledge, the first case is a rare example where we 
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can provide the exact value of the energy gap for an interacting spin chain with alternating 
bond strength. It should be emphasized that this result is both at a feasible value of the 
spin-Peierls dimerization and in a physically meaningful regime of the coupling V The 
utility of these exact results are threefold. First of all, since the lowest gap excitation for case 
(I) is very small, it provides us a relevant check of the reliability of numerical and approximate 
non-perturbative methods to distinguish between gapped and gapless phases. Moreover, we 
present evidence that the second case is an example of a system with a spontaneous spin-Peierls 
effect. This lends support to the prediction by Nersesyan and Tsvelick that negative four- 
spin interaction can induce non-Haldane spin-liquid behaviour in spin ladders. The novelty 
here is that this fact is established at the very strong coupling regime V/J= —4, which is 
beyond the reach of the approach of ref. 0. Finally, together, they motivate us to study the 
gap behaviour for arbitrary values of dimerization and interaction and to search for possible 
gapless phases. 

The exact integrability is derived by identifying the lattice statistical model whose row- 
to-row transfer matrix commutes with Hamiltonian (1) with the parameters 6 and V/J given 
by (2). Such classical statistical system consists of two isotropic six vertex models coupled by 
the total energy-energy interaction. The row-to-row transfer matrix T(A) of a vertex model is 
usually written as the trace of an ordered product of Boltzmann weights, 

T(A) = Tr^[4,!4(A)4'li,^(A) ■ ■ ■ ci'},{X)Cf^{X)] (3) 

where £^'^(A) denotes the local Boltzmann weights with auxiliary space A and quantum space 
j = 1, ■ ■ ■ , L, parametrized by the spectral parameter A. The auxiliary and quantum spaces 
correspond to the horizontal and vertical degrees of freedom of the coupled six vertex models. 
In our case there are four possible states per bond. 

Essential to our approach is to notice that the Boltzmann weights can be conveniently 
written in terms of two commuting Temperly-Lieb operators Ej"^ = 2{Sj"\S^^ — Jj_4/4), 
where /^^^ is the identity operator. More precisely, defining £j-^(A) = Pj,ARf^A{^) with Pj^^ 
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denoting the exchange operator, we found that the i?-matrix Rj°']^{X) is given by 

<i(A) = Ij,A + ^i(A)i?S + u;,{X)E^l},Ef^ (4) 

The weight uj2{X) plays the role of the four-spin interaction V while a combination with 
ci;i(A) is responsible for the dimerization 6. These weights belong to a more general class of 



systems, that originally appeared in integrable coupled Potts models |T6[, and were recently 
rediscovered in the context of Lorentz lattice gases and Fuss-Catalan algebras |T^. For the 
cases we are interested in, their explicit expressions are 

1 for (l) — ^ (I) 

^,(A)= 2 ^' and i^j(A) = 3",?;^, (5) 

where 7 = ln(2 + VS). 

As usual, the corresponding Hamiltonian is obtained as the first-order expansion in A of the 
logarithm of T(A). By using the weights (5) it is not difficult to verify that we indeed recover 
the Hamiltonian (1), up to irrelevant re-scaling constants, at the values given by (2). To make 
further progress we have to explore other properties of the transfer-matrix. In particular, we 



are interested to establish the inversion relation [18, 19, 20 1, since it in principle provides us 
the means to compute the ground state energy and excitation properties. We will start by 
considering the most involved model which is the case with non-null dimerization 6 = 1/2. 
For this system, the inversion identity follows from a combination between the usual unitary 
property, jC^i^2W'^2^i{~-^) = 11,2, and a less standard crossing relation for the weights >Ci"2(A). 
These operators are not independently crossing symmetric, instead they satisfy a novel "mixed" 
crossing property, which reads 

/:S5(A) = M [/:&"Y^(ln3 - A) M, a = 1,2 (6) 

k 

where denotes the transpose and M is a 4 x 4 antidiagonal matrix both acting on the space 
k. 
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Together with unitarity, this crossing relation allows us to derive the following inversion 
identity, 

"cui(ln3 + A)'^ 



Id + T+(A) (7) 



where T~^{X) is a matrix whose elements for large L are exponentially small, 0(e~^), compared 
to the term proportional to the identity Id. This means that, in the thermodynamic limit, the 
last term in (7) vanishes, providing us a much simpler functional equation for all the transfer 
matrix eigenvalues. In particular, the largest eigenvalue AW(A) per site satisfies 

A£(A)Ag(ln3 + A) = "^(^"' + ^^ (8) 

With the help of unitarity, A(^j(A)A^^](— A) = 1, it is possible to solve the functional equation 
(8) under plausible analyticity assumption in the region < A < ln3 where the Boltzmann 
weights are positive. The solution is given by 



We now have the basic ingredients to derive the exact value for the ground state energy per 
site E^g) / J of Hamiltonian (1) at 5 = 1/2 and V/ J = —2. This value is obtained by computing 
Ef}/J = -^lllA^l^^Q + 3/8, which reads 

5 °° 3^"' 
4iV^ = -g - 12E (32,>i_i)(32,+2_i) = -1-433126534 (10) 

The inversion relation and trigonometric periodicity impose stringent constraints to the 
form of the low-lying excitations. They should be described in terms of meromorphic functions 
having two independent periods 21n(3) and 2tt. This observation alone enables us to calculate 
the dispersion relation pO[ and from that one obtains the exact value for the gap. Such gap 
corresponds to the energy necessary to create an excitation with total spin S"^ = 1 in the 
spin-orbital model (I). Here we omit further technicalities and present only the final result. 
We found that the triplet energy gap is given by 



32i+i - 1) re-^32J+2 - 11 



A (I) OO 

i=l L 



:i-3-^/^) 

[1 + 3-^/2) 
4 



2 



0.002869614 (11) 



Interesting enough, the energy gap is very small and this has an immediate application. It 
could be used to test if a given numerical or approximate non-perturbative method can really 
make a clear distinction between a small gap and a real gapless phase. In general, this is a 
difficult task, and we expect that our exact result will be quite relevant to determine suitability 
of multiprecision methods in spin ladders models. For instance, this should be important when 



one wants to predict the scaling behaviour of the g function of the dimerization p2i. 

We turn next to the second solvable point 6 = and V/J= —4. In this case, uJi{X) is null 
and we are only left with the product of two commuting isotropic six vertex models. From the 
point of view of the classical statistical model, this means that we are in fact dealing with an 
alternative representation of the sixteen-state Potts model. By now several properties of the 
general g-state Potts are fairly well understood. In particular, the ground state and the triplet 
gap can be determined either by the inversion trick as above |^ or by a direct mapping onto 
the XX Z Heisenberg chain ||2l|. Their exact expressions are 
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J 

and 



^(11) 

--V3 



1 + 4Et- 

.7 = 1 



(2 + V3)2J 



A (II) 

2v^n 



[1 - (2 + v^)- 



-I 2 



-1.984444091 (12) 
4 ^ ^ 



0.779604542 (13) 



.[l + (2 + v^)-J] 

Further interesting results can still be derived from the mapping of the homogeneous spin- 
orbital model (1) at V/J = —4 onto the the antiferromagnetic XX Z chain with anisotropy 
jz I jx _ 2 First, it is possible to show that the gap of the ffist excitation in the sector 

of total spin S*^ = vanishes exponentially as L ^ oo. The momentum of this excitation is 
TT, which is compatible with the interpretation that in the thermodynamic limit the system 
has two spontaneously dimerized ground states, in accordance with the Lieb-Schultz-Mattis 
theorem ||23|, |2^. This result supports the prediction by Nersesyan and Tsvelick |0 that a 
four-spin interaction may induce dimerized phases in spin ladder models. The bosonization 
arguments of ref.|^ for the weak coupling regime << 1 together with our exact result 

at the very strong coupling point V/ J = —4 indicate that such dimerized phase should be 
robust for a rather large region of ^ < 0. Next, one can explain the numerical observations 



by Pati, Singh and Khomskii P] that the model (II), now for J < (point B in figure 1 of 
ref. 0), possesses an infinitely degenerated ground state. In fact, J < corresponds to the 
ferromagnetic regime of the XX Z chain, whose T = finite entropy for J^/ = —2 is exactly 
computed to be ln(2 + v^)- This value is a rigorous confirmation of the lower bound proposed 
in ref. for such residual entropy. 

Considering these exact results, it is natural to ask if a combined effect of dimerization and 
negative four-spin interaction could lead to a gapless regime. To investigate this problem we 
numerically diagonalize the Hamiltonian (1), with periodic boundary conditions, up to L=14 
sites by using a Lanczos-type algorithm. In table 1 we exhibit our numerical results for the 
ground state energy and the triplet energy gap for both cases (I) and (II). The extrapolations 
towards the infinite volume limit was performed by using the Van den Broeck-Schwartz method 
of convergence [^. While the results for the ground state are in good agreement with the exact 
values, the gap estimates have a rather poor accuracy. This emphasizes the importance of our 
exact results, since they clearly show that one cannot trust the numerical gap estimates beyond 
two significant digits. However, this does not prevent us to make a qualitative comparison of 
the behaviour of the gap in the whole dimerization region < S < 1 for various values of 
the interaction V. This is illustrated in figure 1 for three values of V in the strong coupling 
regime. We observe that indeed for each value of V/ J the gap has a very small minimum for an 
appropriate value of the parameter 6. Note that for V/J= —2 this minimum occurs very near 
the solvable point 5=1/2 where the gap is small but still finite. This helps us to establish an 
upper bound for V/ J, beyond which one probably should rule out strictly null mass gaps. For 
V/J < —2 we observe, however, that this minimum decreases faster towards zero, suggesting 
the possibility of a gapless line on the plane {6, V/ J) in the < V/ J < —2 regime. It would be 
interesting to confirm the existence of this massless line via more powerful numerical methods 
such as the density matrix renormalization group ^6). In particular, this method will allow 
us to determine with good accuracy the critical value V/ J where the gapless line starts. In 
practice, however, when a gap is as small as we calculated, its effect would be invisible at even 
low temperatures which may still be considerably higher than the value of the gap. 
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In summary, we have pointed out the existence of two integrable points in the spin-orbital 
model where relevant physical quantities such as the ground state energy and the excitation 
gap can be evaluated exactly. These results together with the solvable SU(4) symmetric point 
seem to exhaust all possible Bethe ansatz integrable cases of Hamiltonian (1). In addition, our 
numerical analysis indicate the possibility of a gapless line on the plane {6, V/ J) for <V/J < 
—2. This then will add other example of coupled spin chains in which suitable combination of 
dimerization and interaction strength is capable to close the energy gap ||. We also hope 
that our observations will motivate further numerical and analytical investigation in related 
systems. 
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L 








A(")/J 


4 


-1.625 


-2.25 


2.267 949 


4. 


6 


-1.510 455 


-2.083 333 


1.529 642 


2.876 894 


8 


-1.474 743 


-2.033 159 


1.162 347 


2.306 000 


10 


-1.459 071 


-2.012 015 


0.939 894 


1.957 187 


12 


-1.450 822 


-2.001 389 


0.790 102 


1.722 250 


14 


-1.445 954 


-1.995 437 


0.682 160 


1.553 984 


Extrap. 


-1.433 1(±1) 


-1.983 (±1) 


0.011 (±3) 


0.744 (±2) 


Exac-t 


-1.133 126 


-1.981 111 


0.002 869 


0.779 601 



Table 1 : Finite size and extrapolated results for the ground state energy and the triplet 
gap for cases (I) and (II). For comparison we also exhibit the exact results. 
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Figure 1: Extrapolated triplet energy gap A/ J for V/J = —1, —2, —3. For F/J < — 4 we no longer observed 
a local minimum for the gap in the interval < 5 < 1. 
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